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Multi-object filtering with stochastic populations 
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Abstract —While the design of automated knowledge-based 
sensor scheduling is relevant to many multi-target detection 
and tracking problems, tracking algorithms are rarely built for 
this purpose and their outputs provide little flexibility for the 
design of sensor policies. In this paper, we present an estimation 
framework for stochastic populations in the context of multi¬ 
target estimation problems. Fully probabilistic in nature, it allows 
for the evaluation of the population of targets through statistical 
moments, as well as the assessment of sensor observations 
through information-theoretical gain functions. We present a 
principled solution derived from this framework addressing 
challenging multi-target scenarios involving missed detections 
and false alarms, the filter for Distinguishable and Indepen¬ 
dent Stochastic Populations, which propagates information on 
previously-detected targets as well as yet-to-be-detected targets 
while maintaining track continuity. 

I. Introduction 

Multi-target detection and tracking problems, often con¬ 
strained by limitations in sensor capabilities in realistic sce¬ 
narios, are a source of inspiration and motivation for the 
design of automated knowledge-based sensor scheduling poli¬ 
cies optimizing the allocation of sensing resources. Initial 
approaches [ 1 ]—[3] to the multi-target problem, including the 
Multiple Hypothesis Tracking (MHT) and Joint Probabilistic 
Density Association (IPDA) techniques [4], were designed as 
extensions of the single-target tracking problem to the multi¬ 
target case. They rely on an intuitive approach and heuristics 
in order to represent the uncertainty on the number of targets 
in the population of interest, including ad-hoc mechanisms for 
track creation and deletion. Examples of exploitations of these 
track-based approaches in sensor scheduling problems can be 
found in [5], [6], or in [7, chap. 14]. The more recent Finite Set 
Statistics (FISST) framework [8], on the other hand, is fully 
probabilistic in nature. It allows for the principled derivation of 
filtering solutions such as the Probability Hypothesis Density 
(PHD) filter [9], or the Labeled Multi-Bernoulli (LMB) [10] 
filter, and paved the way for the construction of principled 
sensor scheduling policies fully integrated to the multi-target 
estimation problem [11]—[15]. 

In this paper, we present an estimation framework for 
stochastic populations [16] as an original approach to multi¬ 
target estimation problems, and a principled filtering solution 
derived from this framework, the filter for Distinguishable 
and Independent Stochastic Populations (DISP). Some of the 
key features of this approach highlighted in this paper are as 
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follows: 

1. Stochastic populations represent uncertainty on the 
composition of the population of targets and on the targets’ 
state. Two types of statistical moments are presented for a 
stochastic population, in order to evaluate the composition of 
the population or the state of its members. The latter extend 
the regional statistics for FISST-based filters [17], exploited 
in [14] for a sensor scheduling problem. 

2. The estimation framework uses a common probabilistic 
description for a subpopulation of targets indistinguishable 
from each other for the purpose of estimation (e.g. the targets 
that appeared at the same time step and have not been 
detected yet), so that unnecessary permutations are avoided 
during data association steps with observations collected from 
the sensor. 

3. The DISP filter propagates information on previously - 
detected targets, for which track continuity is maintained 
through their past observations, and on yet-to-be-detected as 
well, based on the operator’s knowledge about the targets 
entering the surveillance scene. 

4. Principled information-theoretical gain functions are 
derived from the output of the DISP filter, focusing on 
specific targets and/or specific regions of the surveillance 
scene, in order to quantify the information gain provided by 
the observations collected from the sensor system. It extends 
the information gain functions presented in [18]. 

The paper is organized as follows. Section II gathers all the 
notations used throughout the paper. Section III presents the 
estimation framework for stochastic populations, and its salient 
features in the context of multi-target estimation problems. 
Section IV then presents a typical multi-target surveillance 
activity, and the modelling assumptions leading to the DISP 
filter. Section V then shapes the modelling assumptions into 
the stochastic population propagated by the DISP filter, and 
Section VI and Section VII describe the time prediction and 
data update steps of the DISP filter, respectively. Section VIII 
discusses the exploitation of the DISP filter to assess the pop¬ 
ulation of targets, and introduces a principled approximation 
of the filter with reduced complexity. Section IX concludes, 
followed by additional material and proofs in the appendices. 

II. General notations 
A. Measure theory 

We assume that random variables are defined on some 
probability space (fi, J 7 , P) and denote E the expectation with 
respect to (w.r.t.) P. We denote by M(E) (resp. V(E)) the set 
of finite positive measures (resp. probability measures) on a 
given measurable space (E, 13(E)), where B(E) is the Borel 
u -algebra on E. 
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We shall also use the notation da;, defined for any point 
x £ E, for some infinitesimal Borel set in B(E) containing 
x. In the course of the paper, we shall consider augmented 
spaces of the form X = {ip} U X. The notation da; is then 
to be understood as the Borel set {ip} if x = ip, or as some 
infinitesimal Borel set in £>(X) containing a;, if a; £ X. 

Additionally, the Banach space of all bounded and mea¬ 
surable functions on E equipped with the uniform norm will 
be denoted £°°(E), and we write p(f) = J p(dx)f(x) for 
any /a £ A4(E) and any function / £ C°°(E). For any 
[i,\ji! £ M(E), we shall use the notation p(dx) = //(da;) 
when p(f) = //(/) for any function / £ C°°(E). 

Finally, for some subset B £ 13(E), the indicator function 
1 b is the function on E defined by 


III. Stochastic populations for multi-target 

ESTIMATION PROBLEMS 

Here we provide the salient features of the general estima¬ 
tion framework for stochastic populations [16] that we shall 
exploit for a multi-target estimation problem. For the rest 
of the paper, we denote by X a population of individuals 
of interest, or targets, on which some operator wishes to 
determine individual characteristics gathered into a target state 
(position, velocity, etc.), and we shall assume that the targets 
behave independently from each other. In this context, one can 
describe two sources of uncertainties regarding the population 
of targets: its composition, and the individual state of its 
members. 


1 b(x) 


1, x £ B, 
0, x £ B, 


(1) 


and the scalar p(1b) will sometimes be denoted p,(B), both 
notations being used interchangeably in the rest of the paper. 


B. Time prediction 

A transition kernel Q : E x B(E') —» R from a measurable 
space E into a measurable space E' verifies i) i £■ Q(x,B) 
is measurable for any B £ B(E') and ii) B Q(x,B) is a 
measure for any x £ E. If Q(x, E') = 1 for any x £ E, then 
Q is referred to as a Markov kernel from E to E'. In this 
case, let Tq be the following change of probability measure: 

r Q : V(E) -A V(E’) 

T ^ Fq(/x), 

where rQ(/r)(dx') = p(Q(-,dx')). 

If Q denotes the Markov transition kernel from a space 
Et-i representing the target state at time t — 1 to a space Et 
corresponding to time t, then Yq transforms the information 
maintained on the state of some target from time t — 1 to t. In 
a linear Gaussian problem, it describes the usual single-target 
prediction step of the Kalman filter [19]. 


C. Data update 

Let G £ C°°(E), the following change of measure is 
referred to as Boltzmann-Gibbs transformation [20]: 


■$> G : M(E)^V(E) 
r/i-A 


(3) 


where, assuming 17 (G) > 0, 


T' G (? 7 )(da;) 


1 

V(G) 


G(x)r](dx). 


(4) 


If G denotes the likelihood function for some observation z 
produced by a sensor, then T q transforms the information 
maintained on the state of some target at time t once it is 
associated to z. In a linear Gaussian problem, it describes 
the usual single-target/single-observation update step of the 
Kalman filter. 


A. Multi-target configuration measures 

Assume that E denotes a target state space suitable for the 
description of individual target states, and that the operator 
possesses a set {pi : i £ 1} of probability distributions in 
V(E) indexed by some set I, such that pi = pj => i = j 
for any i,j £ I. The track (i,Pi) associates a unique target 
index j £ I to a corresponding description of the target’s state 
in E through the probability distribution (or law ) pi of the 
target. In the context of filtering, as it will be seen with for 
the DISP filter in this paper, the composition of the index set 
I and the corresponding laws are enriched along the scenario, 
most notably when new observations are collected from the 
sensor system observing the scene. 

We define a multiplicity on I as a family of integers n £ N 1 
indexed by I, describing a composition of the population 
X where there are n, individuals represented by the track 
(■ i,pi), i £ I. It is shown in [16, chap. 2] that a suitable 
representation of the corresponding population is given by the 
multi-target configuration measure p n £ A4(V(E)), a integer¬ 
valued measure on the space of probability distributions V(E) 
defined as 

Mn = ^n»<5[Pi], (5) 

ie 1 

where 5\p\ is the Dirac measure at point p. The multi¬ 
target configuration measure p n embodies all the information 
maintained by the operator on the population X, given the 
specific composition n. The case n, > 2 indicates that several 
targets are indexed by the same track (i,Pi): the law pi 
describes collectively the state of each target of the population 
indexed by i, and these targets are indistinguishable for the 
purpose of estimation. Conversely, the case rij = 1 indicates 
that a single target of the population is characterised by the 
track ( i,pi ), and is thus distinguishable for the purpose of 
estimation. The concept of target distinguishability is a key 
element of the estimation framework for stochastic populations 
[16], and plays an important role in the construction of the 
DISP filter. 

While the form (5) is concise and convenient to describe 
the propagation of information in the context of filtering, it 
is shown in [16, chap. 2] that //„, induces an equivalent joint 
probability measure which enables events regarding the state 
of individuals to be assessed. 
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Example 1. Assume I = {a, b,c}, the target laws p a , Pb,Pc, 
and the multiplicity n £ N 1 such that n a = 1, rib = 2, 
n c = 0. The corresponding multi-target configuration (i n is 
then 

dn = S[p a ] + 2S[p h }. (6) 

Consider some subsets A,B £ B(E) as regions of the target 
state space in which one wishes to assess the presence of 
targets. 

Then, according to p n , the probability that there is one 
target with index a lying within A, two targets with index b, 
one lying within A and one within B, and no target with index 
c, is 

Pa.(A)p h (A)p h (B), ( 7 ) 

and the probability that there is one target with index c is 
equal to zero since the composition of the population does not 
agree with the multiplicity n. 

The formal correspondence between a multi-configuration 
measure p rl and the induced joint probability measure is 
provided in Section A of the Appendix, though it is not 
essential for the rest of the paper as the form (5) will be 
exploited from now on. 

B. Stochastic population 

We have seen in Section III-A that the information on the 
population for a given multiplicity n is represented by the 
multi-target configuration measure ji n . The uncertainty on the 
composition of the population is then embedded in a random 
configuration measure £, i.e. a random variable on the set of 
integer-valued measures 1 on 'P(E), characterised by 

C{F) = Y / N i F(p i ), ( 8 ) 

iei 

for any F £ C°°(P(E)), where N is a random variable on N 1 
induced by £. The stochastic population £ thus encapsulates 
all the uncertainty about the population, as N describes its 
composition and the laws pi describe the states of its members. 
Denoting by c the probability mass function of N, the law of 
the stochastic population £ is given by 

P<s. = c{n)6[p n ], ( 9 ) 

nGN 1 

i.e., a realization of the stochastic population £ is the multi¬ 
target configuration //„, given in Eq. (5), with probability 
c(n). 

Remark 1. The case where the number of target represen¬ 
tations is reduced to one (i.e., |I| = 1 ) is interesting to 
consider. The stochastic population ( 8 ) then reduces to a 
random variable N on N, describing the number of targets 
in the whole population X, and a single law p £ V(E) 
describing the state of each of the targets; in other words, 

[n a more general case beyond the scope of this paper [16], one might 
not be able to index the targets through their probability distributions, i.e., 
the countable index set I might not be available. In this case, the stochastic 
population £ is a general point process on V(E) which cannot be reduced to 
a random multiplicity as in Eq. (8). 


the stochastic population £ becomes equivalent to an inde¬ 
pendently identically distributed (i.i.d.) point process on E 
[21]-[23]. 

C. Population and statistical moments 

Similarly to usual random variables, we can produce 
statistical moments of the stochastic population £. Since a 
stochastic population represents two levels of uncertainty 
regarding the composition of the population and the state of 
its members, we can define moments of two different nature: 
the (full) moments provide statistics on the composition of 
the population (e.g. how many targets are represented by 
some probability distribution p), while the collapsed moments 
provide statistics on the state of the targets (e.g. how many 
targets are lying in some region B of the surveillance scene). 

1) Full moments: The first moment measure 

Mg; £ M(V(E)) and the variance Var<r : C°°(P(E)) —> R 
are defined as 

M e (F)=E[£(F)], (10) 

Var c (F) = E[£(F) 2 ] - E[£(F)] 2 , (11) 

for any F <E C°°(P(E)). 

The component Ni in Eq. ( 8 ), i £ I, is a random variable on 
N describing the size of the subpopulation of targets indexed 
by i. Using the law of the population (9), we can write the 
first moment and covariance of these components as 


m £ (i) = E[7Vj] (12a) 

= ^2 c ( n ) n i . ( 12 b) 

raeN 1 

CO v e (i,j) = E[JVjJVj-] - E[JVi]E[JVj] (13a) 

= ^ c(n)rii [rij - m c (j)], (13b) 

neN 1 

such that the first moment M<j and variance Var e - are found 
to be 

M £( F ) = E m e(*) F (ft)> (14 > 

*€l 

Var c (F) = ^ cov e (i,j)F(p i )F(p j ). (15) 

i,jei 


Since the first moment Me and variance Var c - are statistical 
quantities on V{E), i.e., on the space of probabilities measures 
on the target state space F, they are useful to produce statistics 
on the laws of the targets of the population. 

Example 2. The number of targets represented by a probabil¬ 
ity distribution within some subset B £ B(P(E)) has mean 
and variance given by setting F = 1 b, i.e. 

M £ (1 b) = ^m £ (f)l B (pi), (16) 

iei 

Var c (l B ) = ^2 covc{i,j)l B (j>i)lB(Pj)- (17) 

i 

In particular, considering B = {p,} yields mjfi) and 

coVe(i,f), i.e., the mean size and variance of the subpop¬ 
ulation of targets indexed by i £ I. 
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Also, if a metric d is defined on V(E) then a subset of the 
form B p = {p' £ V(E) \ d(p,p') < 5}, where <5 £ [0,oo), 
can be used to estimate the number of targets whose laws are 
close to p in the sense of d. 

2) Collapsed moments: The collapsed first moment 
mj £ M(E) and the variance var^ on C°°(E) are defined 
as 


m £ (/) = M c (x/), (18) 

var c (/) = Var c (x/), (19) 


where the transformation Xf is given by 

Xf ■ V{E) -> E 

P^p(f), 


( 20 ) 


for any function / £ C°°(E). From Eqs (14), (15), the 
collapsed quantities become 

m <t(/) = Y m c(*)^(/)> (21) 

igl 

var c (/) = Y cav z(i,j)Pi{f)Pj(f)- (22) 

i,je i 


Since the collapsed first moment me and variance v;:»x are 
statistical quantities on the target state space E, they are 
useful to produce statistics on the states of the targets of the 
population. 

Assuming that E is the augmented target state space 
X = {0} U X describing the surveillance scene, where 0 is 
the state of targets absent from the scene (see Section IV-A): 


Example 3. The number of targets lying in some subset 
B £ £>(X) has mean and variance given by setting f = 1 B , 
i.e. 


m c (l B ) = Y m e{i)Pi(B), (23) 

»ei 

var c (l B ) = Y cov €{iJ)Pi{B)pj(B). (24) 

i,j 61 

In particular, considering B = {0} yields the mean and 
variance of the number of targets absent from the scene. 


Remark 2. Moments can also be exploited to evaluate events 
regarding the states of specific subpopulations of£. For exam¬ 
ple, the number of targets lying in some region B £ B(X.) and 
belonging to a subpopulation indexed by either i or j, where 
i,j £ I, has mean and variance given by Mc(l{ Pi Pj }.Xi B ) 
and Varc (1 } Xis )> respectively. 

Note that the collapsed mean and variance produced in 
Example 3 are equivalent to the regional statistics for point 
processes, introduced in [17] in the context of FISST-based 
filters. We will see in Section VIII that the statistical moments 
can be exploited in a similar way from the output of the DISP 
filter. 


D. Merging of populations and mixture of laws 

The number of indexed subpopulations in £ reflects the 
complexity of the information maintained by the operator on 
the population X; in the context of filtering with stochastic 
populations, the computational cost of the tracking algorithms 
is likely to increase with the number of subpopulations |I| 
through the data association phase (see Section VII). For this 
reason, it is sometimes of interest to merge two (or more) 
subpopulations, i.e., replace them by a single subpopulation 
in £, whose probabilistic description accounts for those of the 
subpopulations about to be merged. 

Suppose that the operator wishes to merge two subpopula¬ 
tions with respective index i,j £ I, where i ^ j. Denote by 
i(Bj the index of the merged subpopulation to be constructed. 
From Eq. (8) it is straightforward to construct the probability 
mass function c,®^, accounting for the size of the merged 
subpopulation, as 


= P(£(l{ Pi , P3 .}) = n) (25a) 

= P(A0 + Nj = n) (25b) 

= Y (25c) 

71 (EN 1 
rii+rij =n 

for any n £ N. We now need to determine the mixture 
Pi®j £ B{E) of the probability distributions Pi,Pj £ V{E), 
describing the state of the targets in the merged subpopulation 
and accounting for the laws pi and pj. 

Assume that the merged subpopulation has size n, and 
consider the probability distribution pY £ V{E) given by 


(n) / f \ _ I Nj + Nj — n] 

E [£(l {pi , w} ) I N i + N j = n] 

OC Y C ( n ) t(/) + n jPj (/)] > 

nGN 1 
rii -\-Thj =n 


(26a) 

(26b) 


for any / £ C°°(E). The law averages the laws Pi,Pj 
about to be mixed; note in particular that if the number of 
targets originating from one of the subpopulations (say i ) is 
preponderant in the possible compositions of the merged sub¬ 
population of size n, then the influence of p, is preponderant 
in the composition of the mixture p^j. 

In some situations, it is useful to simplify the structure of the 
merged subpopulation i © j by approximating it by a simpler 
i.i.d. representation collapsing the construction (26) over the 
possible sizes, i.e. 


Pi®j(f) 


E [ C ( 1 {poP ; ,'})] 

oc m € (i)pi(f) + m £ {j) Pj (f). 


(27a) 

(27b) 


Merging operations are common in the practical imple¬ 
mentations of target tracking algorithms, and we will see in 
Section VIII that they can be exploited on the DISP filter in 
order to simplify its structure towards an approximate version. 
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IV. Multi-target estimation with the DISP filter: 

GENERAL CONCEPTS AND MODELLING ASSUMPTIONS 

Here we provide a description of a typical surveillance 
scenario for which the DISP filter has been developed, through 
i) the interactions of the targets with the surveillance scene 
(Section IV-A), ii) the nature of the a priori information avail¬ 
able to the operator (Section IV-B), and iii) the interactions of 
the sensor system with the surveillance scene (Section IV-C). 

A. Surx’eillance activity 

In the context of this paper, an operator wishes to gather 
information about some population X of targets (e.g. vehicles) 
while they lie in some defined region of the physical space 
called the surveillance scene (e.g. the surroundings of a 
facility) and over a given period of time. While the existence 
of each individual in X is assumed certain throughout the 
scenario, its presence in the surveillance scene is not: any 
target may possibly enter/leave the surveillance scene during 
the scenario, such that its presence or absence in the scene at 
any time is unknown to the operator and part of the estimation 
problem. 

The time flow is indexed by some integer t. At any time 
t > 0, each target within the surveillance scene is represented 
via a state x belonging to some target state space X t , 
describing physical and measurable characteristics which are 
unknown and of interest to the operator. The nature of the 
target space space depends on the nature of the surveillance 
activity, but is assumed to be a bounded subset of R dt , such 
that any state x £ X t is a real vector of (finite) dimension 
dt, with a mixture of continuous components such as position 
and/or velocity coordinates, and discrete components such as 
a vehicle type, a level of threat, etc. Conversely, the indi¬ 
viduals currently absent from the surveillance scene assume 
the “empty state” t/>; in consequence, each target in X is 
represented by some point in the augmented target state space 

x t = W}ux t 2 

The goal of the surveillance activity is then for the operator 
to estimate, at any time t > 0, the state of every target of the 
population X in the augmented state space, i.e.: i) whether the 
target is currently in the scene, and ii) if so, its state. We shall 
consider the following assumptions: 

Modelling assumptions. Targets in the population X: 

(Ml) behave independently; 

(M2) enter the scene at most once during the scenario; 

(M3) all have state t/j before t = 0; 

(M4) their evolution is Markovian between two successive 
time steps. 

Assumption (M2) implies that if the same target enters twice 
in the surveillance scene during the scenario, it will be treated 
as two different objects for estimation purposes. Following 
these assumptions, we can divide the population at any time 
t > 0 as 

X t = xf UXfUXf, (28) 

2 The nature of the state space may evolve across time if the operator wishes, 
for example, to estimate additional characteristics or modify the physical 
boundaries of the surveillance scene at some point during the scenario. 


where i) Xf are the targets who have not entered the scene 
yet (and may never do so), ii) Xfi are the appearing targets, 
i.e., those who have just entered the scene, iii) Xf are the 
persistent targets, i.e., those who have entered the scene in a 
previous time step (and may have already left it). 

B. Prior information 

At any time t > 0, the operator may possess and exploit 
some information on the population X(' of appearing individ¬ 
uals (e.g. individuals are likely to appear in groups, alongside 
roads or from the North, etc.). We shall assume that: 

Modelling assumptions. Prior information (time t > 0); 

(M5) the prior information on the appearing targets Xf 
can be described by a single population of indistin¬ 
guishable individuals, independent from the persistent 
targets Xf. 

Assumption (M5) states that the targets appearing at a 
common time are indistinguishable; they will remain so until 
detected by the sensor system. Lifting this assumption would 
allow for a more refined description of appearing targets, and 
lead to a more general filtering solution. 

C. Sensor system and observation process 

At any time t > 0, the surveillance scene is observed 
by some sensor system, providing information on the targets 
through a collection of observations (or measurements). An 
observation 2 belongs to some observation space Z t and 
describes physical quantities measurable by the sensor system 
and relevant to a target’s individual characteristics of interest 
(e.g. range, bearing, angle velocity, etc.). The observation 
space is assumed to be some bounded subset of K d ‘, such that 
any observation z £ Z f is a real vector of (finite) dimension 
d' t . The observations may have a mixture of continuous and 
discrete components in the most general form, although a 
discrete observation space is sufficient to interpret the output 
of some sensors (e.g. resolution cells for a radar). 3 The 
measurements collected by the operator may also contain 
spurious observations or false alarms. We shall assume that: 

Modelling assumptions. Observation process (t > 0): 

(M6) a target produces at most one observation (if not, it 
is a missed detection); 

(M7) an observation originates from at most one target (if 
not, it is a false alarm); 

(M8) targets outside the scene produce no observations; 

(M9) observations are produced independently; 

(M10) observations are distinct, and their number is finite; 

(Mil) an obsen’ation produced from a target depends only 
on its current state. 

Assumptions (M6) and (M7) are a central element of the 
DISP filter: they state that individual observations provide 
specific information about individual targets, and that targets 

3 The observation space may not be constant throughout the scenario 
depending on the context, e.g. if sensors of different nature are exploited 
at different times. 
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become distinguishable through sensor detections (see dis¬ 
cussion on target distinguishability in Section III-A). Lifting 
Assumption (M7) means that targets never become distin¬ 
guishable and leads to the construction of a more general 
filtering solution [16, chap. 3], 

Assumption (M10) guarantees that the observations col¬ 
lected by the operator at any time t > 0 can be described by 
a (possibly empty) observation set Zt . In order to account for 
the missed detections, it is convenient to introduce the “empty 
observation” (!) and define, at any time t > 0, i) the augmented 
observation set Z t = {^} U Z t , and ii) the augmented obser¬ 
vation space 7i t = {</>}UZ t . Assumption (M8) implies that no 
information is acquired at time t > 0 oil the population xf 
of targets that have not entered the scene yet. In consequence, 
the DISP filter aims at estimating the remaining population 
X t \ Xf, i.e., the targets that have already entered the scene. 
As an alternative to Eq. (28), we can divide the population at 
any time t > 0 as 

X t = Xf U X° U X* , (29) 

where i) X° are the targets that have already entered the scene, 
but have never been detected so far (also called yet-to-be-de- 
tected targets), ii) X* are the targets that have been detected 
at least once (also called previously-detected targets). 

V. DISP FILTER: TARGET REPRESENTATION 

In this section, we highlight the articulation between the 
structure of the DISP filter and the concepts of multi-target 
configuration measures and stochastic populations presented in 
Section III. In this section, t > 0 represents an arbitrary time 
step relevant to the scenario. As established in Section IV, 
the operator aims to estimate the targets that have entered the 
scene so far, relying on prior information on the appearing 
targets Xf and the collected observations Zp during all the 
past times 0 < t' < t. 

A. Observation paths 

We first introduce the set Y t of observation histories, or 
observation paths up to time t, defined as 

Y t = Z 0 X ... x Z t , (30) 

as well as the set Y t of non-empty observation paths defined 
as Y t = Y t \{4> t }, where cj) t is the empty observation path at 
time t, i.e. the sequence made of t + 1 empty observations ()>. 
Some examples of observation paths are illustrated in Figure 1. 


{(/>,(/>, (/>,(/>) = </> 3 
(0, 0,0,2 3 ) 
(0,Zi,0, 0) 

(0, Zl ,<t>, z 3 ) 

0 ) 


Time t. 


z ! z 2 z 3 

• Zl 



1Z3 


B. Previously-detected targets 

As we have seen in Section IV-C each observation (if not a 
false alarm) provides specific information on one target, and 
thus the previously-detected targets X* are distinguishable for 
the purpose of estimation. A target in X* can then be indexed 
by the pair i = (f a ,y), made of its observation path y £ Yt 
as well as its time of appearance in the scene 0 < t a < t 
(also called “time of arrival”) and can be described by some 
associated probability distribution^ £ V(X- t ), describing the 
current state of the target indexed by i and no one else. 

Note that the probability distribution p\ is defined on the 
augmented state space X t . The scalar p' t (Y t ) denotes the 
probability that the target indexed by i is currently in the scene, 
given that it exists; it is called the probability of presence of 
target i. The restriction of the probability distribution p' f to the 
state space X* is called the spatial distribution of the target i. 

The set of all possible indices for previously-detected 
targets is given by 

r t ={(t',y)\yeY t ,0<t’ <t}. (31) 

It is clear from the construction of the observation paths 
(30) that each target in I* is a valid characterization of one 
individual in X*. However, not every subset of I* forms a 
valid characterization of the individuals X*, since any two in¬ 
dividuals whose observation paths are incompatible, i.e., who 
share a non-empty observation 4 , may not exist simultaneously 
without violating Assumption (M7). The subsets of targets in 
I* made of pairwise-compatible targets are called hypotheses 
and form a set H t maintained by the DISP filter 5 . 

Following Eq. (5) and given some hypothesis H £ H t , 
the case where the previously-detected targets X* are charac¬ 
terised by the indices in H is then described by the multi-target 
configuration measure 

M? = £*[Pi]- (32) 

i£H 


C. Yet-to-be-detected targets 

As we have seen in Section IV-B the targets appearing in the 
scene simultaneously are indistinguishable for the purpose of 
estimation, and remain so until their (eventual) first detection. 
The subpopulation of targets that appeared at time 0 < t' < t 
and have never been detected so far, undistinguishable from 
each other, are represented with a single index ( t',<f > t ), asso¬ 
ciated to some probability distribution pf £ 'P(Xt). In 
other words, a single track describes collectively the state of 
each target in this subpopulation. The notions of probability of 
presence and spatial distribution defined in Section V-B hold 
for indistinguishable targets as well. 

Note that yet-to-be detected targets appearing at different 
time steps t',t" are characterised by distinct indices, and 
described by the distinct probability distributions '^ l> and 


Fig. 1: The observation paths at time t = 3, given a sequence 
of collected measurements. The six possible observation paths 
are listed on the right. 


4 In the situation illustrated in Figure 1, for example, (0,0, 0 , 23 ) and 
( 0 , 21 , 0 , 0 ) are compatible, but ( 0 , 0 , 0 , 23 ) and ( 0 , zi, 0 , 23 ) are not since 
they share measurement 23 . 

5 The proof is given in Section VII-D. 
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p(t respectively. The set of all possible indices for yet- 
to-be-detected targets is given by 

l° = {(f',</> t )|0<f'<f}. (33) 


Following Eq. (5) and given a family of integers 
n £ N( = indexed by times between 0 and t, the case 
where the yet-to-be detected targets Xf are composed of n t > 
individuals that appeared at time t', 0 < f < t, is then 
described by the multi-target configuration measure 


Tt 


o <t'<t 


n v 


Pi 


(*'■</>*) 


(34) 


A. Input 

The current input is the law P t -\ of the form (36), repre¬ 
senting the targets who have entered the scene no later than 
at time t — 1. The previously-detected targets are indexed with 
the set I*_] of the form (31), while the yet-to-be-detected ones 
are indexed with the set I[ J _i of the form (33). 

The case t = 0 deserves special mention. Given Assump¬ 
tion (M3) none of the targets have entered the scene so far, 
the set of hypotheses H_i is thus reduced to the singleton 

H_! = {0}. (37) 


Remark 3. The special case rtf = 1 in Eq. (34) represents the 
situation where the subpopulation of yet-to-be-detected targets 
that appeared at some time 0 < t' < t is reduced to a sin¬ 
gle individual which therefore becomes distinguishable (see 
Section III-A). This has no consequence on the construction of 
the DISP filter and, for the sake of simplicity, we shall refer to 
the yet-to-be-detected targets as indistinguishable in the rest 
of the paper. 

D. Law of the population 

Following Sections V-B and V-C, the set of all possible 
target indices is given by 

(35a) 

= {(t\y)\yeY t ,0<f<t}. (35b) 

The probability mass function on H* x Nf describing the 
composition of the population X is denoted w t : for instance, 
w t (II ■ n) is the probability that the population X is repre¬ 
sented by the configuration (H, n), i.e., the probability that 
the previously-detected targets are those indexed in H £ H t , 
together with n t i £ N yet-to-be-detected targets that appeared 
at time t' , 0 < t' < t. Following Eq. (9) the stochastic 
population £ t maintained by the DISP filter is thus described 
by the law 

Pt= J2 w t{ H ^ n ) s [P-t + / 4 *]> (36) 

ff6H, neN, 

and the data flow of the DISP filter is as depicted in Figure 2. 

Pt Z t 



Fig. 2: Data flow of the DISP filter (time t > 0). Pf is a 
model parameter describing the appearing targets (detailed in 
Section VI) and the observation set Z± is collected from the 
sensor system. 

VI. DISP FILTER: TIME PREDICTION 
In this section we provide a detailed construction of the 
Bayes prediction step of the DISP filter, i.e., we describe how 
the law of the population is updated from the posterior P t _ i 
to the prediction P t \t-i (see Figure 2). In this section, t > 0 
designs an arbitrary time step relevant to the scenario. The 
proof of the DISP prediction can be found in [16, chap. 3]. 


B. Modelling 

Given Assumption (M4) the knowledge of the operator 
regarding the evolution of the targets since time i — 1 is given 
by a Markov kernel from the previous target state 

space X f _i to the current one X t . In particular, the function 
{fi>}) describes the transition for a target to the empty 
state fir, since Assumption (M2) implies that no target may re¬ 
enter the surveillance scene, it holds that 

m t - lyt (fij,{fij}) = 1. (38) 

In target tracking applications it is then customary to rewrite 
mt-i.t through a reduced Markov kernel rh t -i,t from X t „i 
to X t and a probability of survival p s<t on X t _i such that 

{ m t -i,t(x,dx') = p s ,t{x)m t -i,t{x,dx'), x £ X t _i,a:' £ X t , 
m t -i. t (x, {fii}) = 1 x £ X t _i, 

dx') = 0, x' £ X t , 

m t -i, t {fit, {fit}) = 1. 

(39) 

The reduced Markov kernel iht-i,t describes the transition 
for a target within the surveillance scene. In most practical 
problems it accounts for the operator’s knowledge about the 
motion of targets in the physical space - presence of roads 
and/or obstacles, maximum speed per vehicle type, etc. The 
probability of survival p Si t(x) is a scalar that describes how 
likely is a target, with state x at time t — 1, to be still in the 
scene at time t. 

Given Assumption (M5), the prior information on the ap¬ 
pearing targets Xfi is represented by a subpopulation of indis¬ 
tinguishable targets indexed with * a = (t,fi> t _ 1 ). Their state 
is described by some probability distribution p*“ t _ 1 £ "P(X t ), 
and their number by some probability mass function p\ on N, 
both parameters of the DISP filter (see Figure 2). Note that, by 
construction, appearing targets are currently in the surveillance 
scene. It thus holds that 

p^ t _ 1 (X t ) = l, (40) 

that is, an appearing target is present in the scene almost surely. 

The case where there are n > 0 appearing individuals is 
thus represented by the multi-target configuration measure 

^It-i = ( 41 > 

and following Eq. (9) the population of appearing individuals 
is then described by the stochastic population whose law 
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is given by 


Pt = J2pt(n)S 

n> 0 


, 77 . 

Mtlt-l 


(42) 


C. Target prediction 

The information gathered so far on any target i £ I t _i 
is described by a probability distribution p l t _ 1 on the former 
state space X t _i; in the prediction step, it is transferred to the 
current state space X t through the Markov kernel rri t _ \ t (39). 
The resulting probability distribution £ V(X t ) i s found 

to be [24] 

Pt\t-i = (43) 

that is, as detailed in Section II-B, 

Pt\t-M x ) = [ rn t -i,t(x',dx)p l t _ 1 (dx'). (44) 

JX t _ i 

Note that the prediction step (44) applies to all targets, 
regardless of their distinguishability. It is formative to write 
explicitly the evolution of the probability of presence of a 
target in the scene: 

Property 1. Let i £ I t _i. The probability of presence of the 
predicted target i is found to be 

Pt\t-i( X t) =$_i(Xt_i) - f [l-ps^x^pl^dx). (45) 

4x t _ i 


E. Output 

Following Eq. (43) an hypothesis H £ H t _ 1 is now 
described by the multi-target configuration measure 

= (49) 

i£H 


while the subpopulation of yet-to-be detected targets composed 
of rif individuals that appeared at time t', 0 < t' < t, is now 
described by the multi-target configuration measure 




E 

0 <t'<t 


Tit'6 


Ptlt-l 


(50) 


The law of the stochastic population £t\t-i maintained by the 
DISP filter thus becomes 


p t\t~i= E E w t\t-i( H ’ n ) s [^\t-i+p?\ t -i]- ( 5l ) 

n6N t 


VII. DISP FILTER: DATA UPDATE 
In this section we provide a detailed construction of the 
Bayes update step of the DISP filter, i.e., we describe how the 
law of the population is updated from the prediction Pt\t -i to 
the posterior P t (see Figure 2). In this section, t > 0 designs an 
arbitrary time step relevant to the scenario, and Z t an arbitrary 
observation set collected from the sensor system. The proof 
of the DISP update can be found in [16, chap. 3], 


The proof is given in Appendix B-A. Given Assump¬ 
tion (M2) no target can re-enter the scene and thus, as 
confirmed by Eq. (45), the probability of presence of a target i 
is non-increasing during the prediction step. It is monotonous 
if and only if (iff) the probability of survival is one over the 
support of p\_i - that is, iff the target i stays in the scene 
almost surely - and it drops to zero iff the probability of 
survival is zero over the support of p , t _ 1 - that is, iff the 
target i leaves the scene almost surely. 


D. Population prediction 

The appearing targets are, by construction, yet-to-be- 
detected; the index set I°_ x is thus augmented as follows 


\\t-i ~ 14—1 U {(^ 4>t-i)}- (46) 

On the other hand, since the observation set Z t is not available 
yet (see Figure 2), neither the distinguishable targets indexed 
in I*_ 1 nor the composition of the hypotheses in H t _i are 
modified by the prediction step. The predicted set of target 
indices is thus 

= I*_i U (47) 


Since the subpopulations of appearing and persistent targets 
are assumed independent, the probability mass function on 
H t _! X N f , for any hypothesis H £ H t _i and any family of 
integers h £ N t , becomes 


w t | 


W(_i(7T, n)pt(n t ), t>0 
Po(h 0 ), t = 0, 


(48) 


where n £ N f _i is the family of integers such that n t ' = hf 
for all 0 <t'<t. 


A. Modelling 

Given Assumption (Mil) the knowledge of the operator 
about the observation process, and in particular the imperfec¬ 
tions of the sensor (measurement noise, false alarms, missed 
detections) is described by i) a non-negative real-valued func¬ 
tion gt{-,z) £ C°°(X t ), interpreted as a likelihood and given 
for any z £ Z t , and ii) a probability of false alarm pf at on 
Z f . In particular, the function gt{f>, ■) describes the observation 
from an target absent from the scene; since Assumption (M8) 
implies that no target can be detected unless it is in the scene, 
it holds that 

\9t(ip,z)=0, z£ Z f , 

= i- 

For any observation z £ Z f , it is then customary to write 
gt(-,z) through a restricted likelihood € £°°(X t ) and 

a probability of detection on X t such that 

{ . 9 t{x , z) = p d: t(x)£ t (x, z), x £ X t , z £ Z t , 
g t {x, f) = 1- Pd,t(z), x £ X t , 

9 t(t/j,z) = 0, z£ Z t , 

9t(i’,4>) = !• 


The restricted likelihood describes the accuracy of the sensor - 
£t(x, z) denotes the likelihood that an observation z originates 
from a target with state x. The probability of detection pd,t{%) 
is a scalar that describes how likely a target with state x is to 
be detected by the sensor 6 at the current time t. 


6 The support of the function p^ t , i.e., the subset of the target state space 
Xt where the probability of detection is non-zero, is called the (current) field 
of view of the sensor. Note that the field of view is included in, but does not 
necessarily equal to, the surveillance scene: if the sensor coverage is limited, 
there are “blind zones” in which targets cannot be observed. 






9 


In the context of this paper, the false alarms are spurious 
observations stemming from the sensor system itself, indepen¬ 
dently of the population X (e.g. malfunctions from the sensor 
device, transmission errors, etc.)- For any collected observation 
z £ Z t , the scalar pf a ,t(z) is the probability that z is a spurious 
observation; conversely, the scalar 1— Pi a ,t{z) is the probability 
that z originates from a target in X. 


The proof is given in Appendix B-B. Given Assump¬ 
tion (M8) only targets in the scene can be detected and thus, 
as confirmed by Eq. (59), a target detected this time step lies 
in the scene almost surely. On the other hand, there is no fresh 
evidence on the presence of targets that are not detected this 
time step, and their probability of presence is non-increasing. 


B. Target update 

Since a new observation set Z t is available, the targets’ 
observation path can be updated with a new data association, 
matching a predicted target to an observation in Z t , and 
their probability distribution can be updated accordingly. By 
construction, the yet-to-be-detected targets are now those who 
remain to be detected following the current observation, i.e. 7 


= (54a) 

= {(t',<p t ) I (t', 0 t _ x ) (54b) 

Conversely, the previously-detected targets are now those 
who have been at least detected once following the current 
observation, whether they had already been detected in the 
past (and hence were already distinguishable), or have just 
been detected for the first time (and hence have just become 
distinguishable), i.e. 

I* = I e i t \t-i,z e z u y.z £ <f t }. (55) 

The updated set of target indices is thus 


= I* U1° (56a) 

= {(t',y:z) | (t',y) £ l t \t-\,z € Z t }. (56b) 

For any target i = ( t',y ) £ Iqt-i, whether it is associated 
to an observation ( z £ Z t ) or a missed detection (z = </>), 
we shall exploit the slight abuse of notation i:z to denote 
the updated target ( t',y:z ) £ It. The probability distribution 
P\\ t -i updates following Bayes’ rule through the Boltzmann- 
Gibbs transformation [20] 

pi :z = ** ( .,*)(pi| t -i), (57) 


that is, as detailed in Section II-C, 


Pl :z ( dx) 


9 t(x,z)p\ | t i(dtc) 


(58) 


It is formative to write explicitly the evolution of the proba¬ 
bility of presence of an updated target: 


Property 2. Let i = (t',y) £ Iqt-i and z £ Z t . The proba¬ 
bility of presence of the updated target i:z = ( t',y:z ) £ It is 
found to be 


z £ Zt, 


pr™ = 


Pt| t iM-,</>)ix t 


i -Ptb-i( x t) +p\\ t -\(9t{-,4>)^ t y 


z = 


(59) 


7 “:” is the concatenation operator on sequences, i.e., (ei,..., e n ) :e = 
(ei,..., e n ,e). 


C. Population update 

Following the construction of the updated targets the com¬ 
position of the population is reassessed, i.e., i) the set of 
hypotheses H t _i is updated to reflect the combinations of 
compatible previously-detected targets, ii) the probability mass 
function w t | t _! is updated to assess the composition of the 
population X. The core of the update step consists in the 
data association, where potential sources of observations are 
matched with the new observation set Z t , and every resulting 
association is assessed. The potential sources of observations 
are: 

• the previously-detected targets X°, indexed by I*_ 1 ; 

• the yet-to-be-detected targets Xf, indexed by I^uj; 

• the clutter, generating the false alarms. 

For the purpose of data association, a clutter generator is mod¬ 
elled for each collected observation z £ Z t \ a s discussed in 
Section VII-A, it produces the observation z with probability 
Pfa,t(z) (in which case z is a false alarm), or it produces 
no observation with probability 1 — pf a ,t(z) (in which case 
z originates from a target). 

1) Data association: Let us fix a pair ( H , n) £ H t _i x N t 
corresponding to a given configuration of the population (see 
Eq. (51)). A data association h then associates the sources of 
observations - the previously-detected targets described by H, 
a number of yet-to-be-detected targets described by n, and the 
\Z t \ clutter generators - to the collected observations z £ Z t 
and the empty observation cj> (see Figure 3). 

Given Assumptions (M6), (M7), the set of data associ¬ 
ations Admz t (H, n) can be defined as the set of tuples 
Zd, Zm, u) such that i) H<\ C H contains the previous¬ 
ly-detected targets that are detected this time step, ii) Z,\ C Z t 
contains the observations associated to these detected targets, 
iii) v is a bijective function from 1 f \ to Z ( \ associating detected 
targets with observations, iv) Z{& = {^fd,t'}*' = o is a family 
of disjoint subsets of Z t \ Z&, satisfying |^fd,t'| < n t > for 
any 0 < t' < t, where are the observations produced 

by yet-to-be-detected targets that appeared at time t' and are 
about to be detected for the first time. 

Note that there is only one way to associate the subset 
of observations Z{A,t' to |^fd,t'| inviduals among the n t / 
composing the subpopulation of yet-to-be-detected indexed by 
(f', (ft- 1), because the latter is composed of indistinguishable 
targets for the purpose of estimation. 

Each triplet a = (H, n, h), where h £ Adm z t {H,n), 
then represents a specific association scheme described by the 
scalar 


Ut = V d,t X «md,t X Ufa,*. 


(60) 
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Previously-detected targets 

(h G H t _i) Yet-to-be-detected targets (to G N t ) Clutter generators ( Z t ) 



Collected observations ( Z t ) Empty observations 


Fig. 3: Data association, for a given configuration of the population (FF, n) and a given set of collected observations Z t . 


where 


«d,*=n n n 


ieH d 


u md,t 


IT ^ (&(■></>))) * 1 td x^[ pi 


0 <t'<t 


i&H\H d 


do not violate the “at-most-one-measurement-per-target” and 
“at-most-one-target-per-measurement” rules given by Assump¬ 
tions (M6) and (M7). One may wonder whether the compo¬ 
sition of the population of previously-detected targets X* is 
represented by H t , i.e., whether the possible configurations 
for X* are the hypotheses H t . 


<,t = n ( 1 -Wa,t(^)) x n Pfa.tO), 
z£Zt\Zf a z€zZf a 

( 61 ) 

where Z la = Z t \ ( Z d U (Jo<t'<t Z M.t') is the subset of 
observations associated with clutter generators, i.e., the false 
alarms (see Figure 3). 

2) Hypothesis update: A given association scheme 
a = (FF, n, h) leads to the construction of a unique updated 
hypothesis FF C I* of the form 

H = S d (H d , v) U S md (FF \ H d ) U Sf d (Z fd ), (62) 

with 


Theorem 1. Define the compatibility relation as the symmetric 
binary relation ~ on the set of obserx’ation paths Y t given by 

(Vy,y G Yt) y= = (z' 0 ,..., z' t ), 

y~y'&[ [zt 1 = z' t ’\ => z t ' = <t>,0 <t' <t\, ( 65 ) 

and the set of pairwise-compatible subsets (or consistent 
subsets) of distinguishable targets as 

Const (I*) = {/ C I* | (Vi, i' G I) 

i = (± i' => y' ~ y"|. (66) 


S d (H d ,v) = {(t',y:v(i)) \ i = ( t',y ) G H d }, 
S md (H \ H d ) = {(t 1 , y-4) | (t',y) G H \ H d }, 

^(^td) = C {(i / , <Pt-i :z ) '■ z G Zf djt i}. 

o <t'<t 

Also, define ft G N t as n t / = n t / — |Zf di t/|, 0 < t' < t, such 
that (FT, h) is the updated configuration under the association 
scheme a = (FF, n, h). Using Bayes’ rule, the probability of 
existence of the updated population is found to be 

w t (FF, n) oc wn)u t a . (63) 

The updated hypotheses H t are constructed by considering 
all the possible admissible association schemes a = (FF, n, h), 
where (FF, n) G H f _i x N f , h G Adm z t {H,n). Following 
Eq. (62), the set of updated hypotheses H f can be written as 

H t - jSd(FF d , v) U S md (H \ FF d ) U S {d (Z id ) \ 

(FF d , Z d , Zf d , v) G Adm Zt (FF,n),FF G H t _i,n G N f |. 

(64) 


Then it holds that 


H t = Const (I*), (67) 

that is, the consistent subsets of previously-detected targets are 
the hypotheses. 

The proof is given in Appendix B-C. 


E. Output 

Following Eq. (57) an hypothesis FF G H f is now described 
by the multi-target configuration measure 

d? = J2 S[pI], ( 68 ) 

i£H 


while the subpopulation of yet-to-be detected targets composed 
of n t ' individuals that appeared at time t !, 0 < t' < t, is now 
described by the multi-target configuration measure 






(69) 


D. Hypotheses H 4 and previously-detected targets I* 

We have seen in Section VII-C that the hypotheses in H, 
aims at representing the subsets of previously-detected targets 
that are pairwise-compatible, i.e., whose observation paths 


The law of the stochastic population <£ t maintained by the 
DISP filter thus becomes 

Pt= w t( H ’ n )6[p? + p?]. (70) 

ff£H t nEN, 
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VIII. DISP FILTER: EXPLOITATION 

The law P t of the form (36) propagated by the DISP filter 
(see Figure 2) describes the stochastic population (£ t defined 
as 

C t (F) = J2NlF(pl), (71) 

»€ It 

for any F £ £°°(7 :> (Xt)), where the probability mass function 
c t of the random variable TV, on N 11 ' is such that 

fn°), J* G H t and n* = l,i €/•, 
c, n = < (72) 

1 ^ 0 , otherwise, 

where /* = supp(n) n I* and n° £ N, is given by 
n£, = n ( t /.0 ), 0 < t' < t. That is, a possible configuration 
n £ N 1 * of the population X has zero mass unless the 
previously-detected targets with non-zero multiplicity i) form 
an hypothesis (they are pairwise-compatible) ii) have a multi¬ 
plicity of one (they are distinguishable). 

As we have seen in Section III, several tools can be 
designed in order to exploit a stochastic population, such as 
the extraction of statistical moments or merging operations. 
We shall see in this section how these tools can be exploited 
on the output of the DISP filter for meaningful operations in 
a multi-target tracking scenario. 

A. Maximum a Posteriori 

The operator may wish to produce a Maximum a Posteriori 
(MAP) estimate of the population X in order to display the 
most probable multi-target configuration from the output of the 
DISP filter. From the law (36) it is straightforward to extract 
the most probable configuration in H, x N t : 

(if, n)* = arg max w t (H, n), (73) 

(H,n) 

that is, the previously-detected targets are most likely to be 
those in H*, and the number of yet-to-be-detected targets that 
appeared at time t' is most likely to be n*, , 0 <t'< t. 

Some individual targets i £ II* might be tentative tracks 
that a “cautious” operator may not wish to display. The 
probability of existence of distinguishable targets, defined in 
Section VIII-B, is a simple statistics assessing the likelihood of 
individual tracks; a simple threshold can be set on the targets’ 
probability of existence to make sure that only the reliable 
ones, with a degree of confidence set by the operator, are to 
be displayed. 

The extraction of the MAP state of each target to be 
displayed is obviously dependent on the chosen implemen¬ 
tation technique for spatial distributions on the state space, 
and is left out of the scope of this paper. More details on 
Gaussian Mixture (GM) or Sequential Monte Carlo (SMC) 
implementations can be found in [25] or [26], respectively. 


B. Statistical moments 

All the quantities related to the full (14), (15) or collapsed 
(21), (22) moments of a stochastic population are produced 
from the statistical moments (12b), (13b) of the indexed 


subpopulations. In the case of the DISP filter, they are found 
to be 


f E E i&I ti 

He H, nGN, 

m f (*) = < H3i (74) 

E E f£l t °, 

IfleH, nCNf 


and 

co v t (i,j) 


E E w *( ff ’ n ) 

He H, n6N, 




= 


E E w * i ' H - n ) ~ m *(■?)]> 

HeUt nGN t 
H3i 


E E n ) n i\n.j - m t(j)], 

, HeHt ngN, 


i,3 G I?, 

* G I*, j £ 1° 
hj £ H° ■ 

(75) 


Note in particular that the expected size m,(T) of a population 
associated to a previously-detected target i £ I*, as shown in 
Eq. (74), is a scalar between 0 and 1: it is equal to i) 0 if 
the target i does not belong to any hypothesis with non-zero 
mass, i.e., if the target i exists almost never, and ii) 1 if the 
target i belongs to all the hypotheses with non-zero mass, 
i.e., if the target i exists almost surely. For this reason, the 
quantity m f (i) is also called the probability of existence^ of 
the previously-detected target i £ I*, and can be exploited to 
assess its credibility (e.g. for displaying purposes, as discussed 
in Section VIII-A). 

Once the statistics (74), (75) of the subpopulations com¬ 
posing the stochastic population £, have been computed, it 
is straightforward to produce meaningful statistics out of the 
full or collapsed moments, as illustrated in Section III through 
Examples 2 or 3. Following Example 3, the regional statistics 
[17] of the population, i.e., the mean and variance of the 
number of targets in any subset B £ T3(X t ), are given by 

= E m 'W p ‘( B )' (76) 

t€ It 

var t (B)= E c °Vt(iJ)p l t(B)p 3 t (B), (77) 

i,je it 

and provide the operator with some information on the level 
of target activity, with associated uncertainty, in the region B 
(e.g. a region of strategic importance in the surveillance scene). 
They can be exploited, for example, in a sensor scheduling 
policy exploring the regions where the uncertainty is the 
highest [14]. 


C. Information gain 

In many applications, notably to establish a hierarchy among 
possible actions in the construction of a sensor scheduling 
policy, some performance metric is necessary in order to assess 


8 Not to be confused with the probability of presence p[(Xt) defined in 
Section V-B, assessing the probability that the target i is still in the scene, 
provided that it exists. 
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quantitatively a sensor action leading to the collection of some 
observation set Z t . We proposed in [18] an information gain 
for stochastic populations, that we shall extend in this section 
to the DISP filter presented in this paper. 

We aim at quantifying the information gain from the pre¬ 
dicted law Pt\t-i of the form (51) to the posterior law P t of the 
form (70) maintained by the DISP filter, given an observation 
set Z t collected by the sensor system. We shall assume that 
the probability distributions p £ 'P(Xt) considered in this 
section admit a Radon-Nikodym derivative with respect to 
some reference measure p, £ A4(Xt), with p({ip}) = 1- For 
the sake of simplicity, the resulting probability density ^ will 
be denoted p. 

1) Target gain: Assume some predicted track indexed by 

1 £ i = some (possibly empty) observation 

2 £ Z t , and the corresponding updated track indexed by 
i:z = ( t',y:z ) £ I* (see Section VII-B). For some mapping 
F : 'P(Xt) —> 'P(Xt), to be specified later on, we define 
the information gain from a target indexed by i, updated with 
observation z, as the Renyi divergence [27] 


om odog / [p(pq t _ 1 )(x)]“ [ F (pi [z )(x)] 


L JX, 


V(dx)J, 

(78) 


where the constant factor is (a — l) -1 , and 0 < a < 1 is the 
order of the divergence 9 . 

2) Configuration gain: Consider now some configura¬ 
tion (if, n) £ H t _i x N t , and some association scheme 
a = (if, n, h), where h £ Adm^ t (if, n) (see Section VII-C). 
Under the association scheme a the configuration (if, n) 
updates to (if, h) £ H t x N t , where the updated hypothesis 
if is described by Eq. (62). We can then define the information 
gain of the association scheme a as 


= E [ E G k*,.,>( F > + K'- i z >«'O 1 

0<t' <t Z £d, t / 

+ ]T Gf\ F )+ J2 G t(F), G9) 

i&H d i<=H\H d 

that is, as the sum of the individual gains (78) from the 
targets in the configuration (if, n). The information gain for 
the configuration follows as 

<W)(F)= £ W aG (Wn ( 80 ) 

hGAdmz ( ( H,n) 


4) Evaluation of the population gain: We can now shape 
the mapping F in order to focus the evaluation of the informa¬ 
tion gain (81) towards specific targets and/or specific regions 
of the state space of interest to the operator [18]. 

Consider the probability distribution £ V(X t ) defined 
as p,/.({/0}) = 1, i.e., all the mass is concentrated on the empty 
target state fi>. For any subset /ci tN , consider the mapping 
F 1 defined as 


F 1 ^- 1) 

F\rt z ) 


Pt\t-n 

i £ I 

(82) 

Pip, 

otherwise, 

P r, 

i £ /, 

(83) 


otherwise. 


The individual gain (78) then becomes 


GKF 1 ) 


OC 


log 


V( dx) 


L «, 


i £ I 

otherwise. 

(84) 


That is, the scalar G^F 1 ) quantifies the information gain, 
provided by the collected observations Z t , regarding the 
targets % £ I. In particular, the scalar Gt(P <: ^ 1 - 1 ) quantifies 
the information gain regarding the whole population. 

For any subset B C £>(X f ), consider now the mapping 
fs : X t —>■ Xt defined as 


Ib(x) 


x, x £ B 

otherwise, 


(85) 


and define Fb(p) as the image of the measure p £ VfX. t ) 
under f b (see Section 3.6. in [28]), or pushforward measure, 
defined by 

(PB(p))(dx) =p(ff, 1 (dx)). ( 86 ) 


The individual gain (78) then becomes 


Gi(F B ) oclog / [p\u 1 (x)} a [pf- z {x)\ la p(dx) 


+ [pi {t _ 1 (X t \B)] a [pf*(X t \B)\ i 


(87) 


where the association weights vf are given by Eq. (60). 

3) Population gain: The expected information gain for the 
stochastic population with law P t \t~i of the form (51), given 
the collected observations Z t , is thus given by 

G t (F) = E[G (H ,„ } (F)] (81a) 

= E E n)G(H,n)(F), (81b) 

KH,_! n 6 N, 

where the expectation in (81a) is taken w.r.t. the configurations 
in H t _i x N t . 

9 The order of the divergence may be dependent on the hack i, but for the 
sake of simplicity we will drop this dependency. 


We see that the gain is zero iff the (possibly empty) ob¬ 
servation z carried no additional information on the target i 
regarding: 

• its localization in B, since p l t ^ t _ 1 = pf z on B, and 

• its presence in B, since pJ| t _ 1 (X t \B) = p l t :z (X t \ B). 

That is, the scalar G^Fb) quantifies the information gain, 
provided by the collected observations Z t , within the region 
B. In particular, the scalar Gt(Fx t ) quantifies the information 
gain in the whole state space. 

In addition, the scalar Gt(Fs o F 1 ) quantifies the informa¬ 
tion gain, provided by the collected observations Z t , regarding 
the targets i £ I and within the region B. 
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5) Knowledge-based sensor policy: The population gain 
(81) extends the results in [18], exploited for an approximate 
version of the DISP filter (briefly discussed in Section VIII-D), 
to the general version of the DISP filter introduced in this 
paper. A similar approach to [18] could then be applied to 
produce a closed-loop sensor scheduling policy focusing on 
sensing actions yielding the highest expected information gain, 
though this is left out of the scope of this paper. 

Note that considering subsets of the form 1 C in the 

expression of the population gain (84), would favour sensor 
actions exploiting previously-detected targets; conversely, con¬ 
sidering subsets of the form I C I° t _ | would favour sensor 
actions exploring the state space for the search of yet-to-be- 
detected targets. 


D. Merging operations 

Merging operations are important for practical implemen¬ 
tations of multi-target tracking algorithms, for they allow to 
mitigate the loss of information whenever approximated ver¬ 
sions of the algorithms are necessary to curtail computational 
costs. We shall highlight in this section a specific merging 
operation leading to a principled approximation of the DISP 
filter whose complexity is significantly reduced 10 . 

In its most general form (36), the DISP filter maintains 
specific representations for populations of yet-to-be-detected 
targets that appeared at different times. If the prior infor¬ 
mation on appearing targets (42) varies little across time, 
however, the probability distributions of the undetected targets 
{Pt\t-i I * ^ It_i} are likely to be close to each other, and 
close to the probability distribution describing the appearing 
targets P^ t _ 1 - In this context, one may wish to merge all the 
subpopulations of yet-to-be-detected targets indexed in I(_ l5 
alongside the population of appearing targets indexed by i a , 
into a single population describing all the yet-to-be-detected 
targets so far, i.e., regardless of their time of arrival. 

Suppose that these merging operations have been conducted 
until the present time t > 0. Since the appearing targets were 
merged in past times regardless of their time of arrival, the 
targets indexed in are characterized with their observation 
path, i.e. 

I?_i = Yt_r, (88) 

I t °_! = {0 t _r}. (89) 


The multi-target configuration (34) for yet-to-be-detected tar¬ 
gets at time t — 1 then takes the simpler form 


/4-1 = nS 


Pt-i 


(90) 


such that the posterior information (36) at time t — 1 takes the 
simpler form 

Pt- 1= +/4-i]- ( 91 ) 

I?6H m nSN 


Using Eq. (25), the merging of the yet-to-be-detected targets 
indexed by 4> t _ l and the appearing targets indexed by i & 


10 Another principled approximation, based on additional modelling as¬ 
sumptions and leading to a much simpler structure, yields the filter for 
Hypothesised and Independent Stochastic Populations (HISP) [16, chap. 4], 


transforms the prediction of the probability mass function (48) 
to the simpler form 




Y w t-i(H,n')pt(n-n), 

0 <n'<n 

Po(n), 


t > 0 

(92) 

t = 0, 


for any H £ H t _i and any n £ N. 

Once the probability distributions {p l t _ 1 \i £ It-i} have 
been transformed to their predicted values {ph t _ 1 | i £ It_i} 

according to Eq. (43), the laws pf* t Z_\ and pf t _ [ can be mixed 
according to Eq. (27). The predicted law (51) then takes the 
simpler form 


Pt\t-i= Y Y w t\t-i{H,n)6[pf ]t _ 1 (93) 

H t _i n£N 


The main simplification occurs in the data association of 
the update step (see Section VII-C). For a given configuration 
(. H,n ) £ H 4 _! x N, the set of data associations Adm.z t (H, n ) 
can be defined as the set of tuples (H^, Z <j, Z^, v) such that 
i) P<\ C H contains the previously-detected targets that are 
detected this time step, ii) Z,\ C Z t contains the observations 
associated to these detected targets, iii) v is a bijective function 
from Hd to Z,\ associating detected targets with observations, 
iv) Zid is a subset of Zt \ Zd satisfying |2Tfd| < n, containing 
the observations produced by the yet-to-be-detected targets 
that are about to be detected for the first time. 

The subset Z{d of observations allocated to first detections 
does not have to be partitioned among the subpopulations 
of yet-to-be-detected targets appearing at different times, re¬ 
sulting in a simplified data association mechanism (see the 
simplified structure illustrated in Figure 4, compared to the 
general structure illustrated in Figure 3). The number of 
hypotheses created in the update step is thus reduced, and 
the number of subpopulations of yet-to-be-detected targets 
in the stochastic population remains fixed at one, rather 
than growing of one at every time step. However, as seen in 
Eqs (88) and (89), the information about the time of arrival 
of the targets has been lost in the merging operation and this 
approximated DISP filter cannot estimate the period of time 
for which a given target has been in the surveillance scene 
prior to its first detection. 

A further approximated version of the DISP filter has 
been proposed, in which appearing targets are assumed to be 
detected at their time of arrival [29]. In this case, since targets 
are assumed detected upon entering the scene, no information 
is propagated across time on yet-to-be-detected targets; in 
other words, this approximated version of the DISP loses the 
ability to estimate the number and states of the targets that have 
not been detected yet. This approximation has been recently 
exploited for a multi-target tracking problem in the context of 
space situational awareness [30]. 
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Previously-detected targets 

Yet-to-be-detected targets (n £ N) Clutter generators ( Z t ) 



Collected observations (Z t ) Empty observations 


Fig. 4: Data association (approximated version), for a given configuration of the population (H , n) and a given set of collected 
observations Z t . 


IX. Conclusion 

In this paper, we present an estimation framework for 
stochastic populations as an original approach for multi-target 
estimation problems, and the derivation of a detection and 
tracking algorithm adapted to challenging scenarios, the DISP 
filter, maintaining information on the targets that have been 
previously detected by the sensor system observing the scene, 
as well as those that have not been detected yet. We show 
that the fully probabilistic nature of the framework can be 
exploited to produce meaningful quantities describing the 
output of the DISP filter, including statistical moments on the 
population of targets or information-theoretical gain functions 
assessing the utility of a set of observations collected from 
the sensor system. This paper focussed on the methodological 
developments of multi-object filters; forthcoming works will 
focus on the illustration of the novel concepts introduced in 
this paper. 


Appendix A 

Multi-target configuration measures and joint 

PROBABILITY MEASURES 

It can be shown [16, chap. 2] than any multi-target con¬ 
figuration measure of the form (5) induces an equivalent joint 
probability measure defined as follows. First, we introduce the 
set E x as 

E x = IJ E k , (94) 

fc>l 

and the set Nj as 

Ni = {(/,n)|/CI,neN J }. (95) 

For a given pair (I, ra) £ Ni, we then introduce the sets 

E$ n) = {X £ (. E x Y | (Vi £ /) Xi £ E ni }, (96) 
E?= IJ E ( r n) . (97) 

(i»eNi 

Then, for a given multiplicity n £ N 1 and denoting by I = 
supp(ro) its support, the multi-configuration //„, induces the 
joint probability measure P(-\I,n) £ V(E X ) defined as 

(VB £ S(4 n) )) P(B\I,n) = Y[p xn *(Bi). (98) 

i£l 


where, for any n £ N and any p £ V(E), p xn £ V(E n ) is 
defined as 

n 

(\/B = (B 1 ,...,B n )&B(E n )) p xn (B) = l[p(B l ). (99) 

2=1 

Note that all the mass of the joint probability measure 
P(jl, to) is concentrated in Ej, i.e, in the elements of E x 
agreeing with the multiplicity to. 

Example 4. Assume I = {a, b, c}, the target laws p a , Pb,Pc, 
the subset I = {a, b}, and the multiplicity n £ N 1 such that 
n a = 1, n b = 2, n c = 0. The multi-target configuration p n 
is then 

hn = <5[p a ] + 2S[Pb\. (100) 

Consider the elements B. 13' £ B(E x ) with B ; , = A, B^ = 
(■ A,B ), B’ a = A, B' c = B. The probability that there is one 
target indexed with a lying within A, two targets indexed with 
b, one lying within A and one within B, and no target indexed 
with c, is then 

P(B\I, n) = p a (A)p h (A)p h (B), (101) 

and the probability that there is one target indexed with a 
lying within A, one target indexed with c lying within B, and 
no target indexed with b, is 

P(B'\I, n) = 0, (102) 

since B’ £ B(E^ n) ). 


Appendix B 
Proofs 


A. Property 1 

Proof: Let i £ I t _i, and let / £ £°°(X t ) be an arbitrary 
function. Substituting the expression of the transformation (2) 
into the definition of the predicted distribution (43) yields 


f(x')m t -i,t(x,dx') 


(103a) 

Pt- i(dz). 

(103b) 


The probability of presence of the predicted track is then given 
by selecting / = lx t in Eq. (103b), i.e. 


Pq t -i(X t ) 






Pt- i(dz)- 


(104) 
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Now, considering the expression of the Markov transition 
kernel (39), Eq. (104) simplifies to 


i4-t(x t ) 


Pe,t(x)lflt-l,t{x, dx') 


/X t 


Ps,t{ X ) 


/ rh t -i,t(x,dx') 

•'X, 


pI-Mx) 


pI-Mx ) 


=1 


(105a) 

(105b) 


= Pt_i(X t _i)+ / [p s , t (a:) - l]Pt_i(dx). (105c) 


B. Property 2 

Proof: Let i = (f',y) e 2 £ Z t , i:z = 

( t',y:z) £ I*, and let / G £°°(X t ) be an arbitrary function. 
The probability distribution of the predicted track i can be 
decomposed as follows: 

Pi\t-i(f) =Pt l t- 1 (lw/)+Pt|t- 1 (lx t /). d06) 

Now, substituting the expression of the Boltzmann-Gibbs 
transformation (4) into the definition of the posterior distri¬ 
bution (57) yields 


C. Theorem 1 

Proof: The proof follows an inductive approach. By 
definition of the original time, no observation were produced 
before t = 0; therefore, the initial set of observation paths 
Y_i is empty. It follows from the definition of consistent 
populations (66) that Const(/* 1 ) is reduced to the “empty” 
hypothesis, i.e., Const(/* 1 ) = {0}. Since, by assumption, the 
set of hypotheses is initialized with the “empty” hypothesis, 
it follows that the base case 

H_i = Const (/*!) (Ill) 

is true. Let us now suppose that the case is true at some rank 
t — 1 > —1, i.e., H t _i = Const(I*_ 1 ), and let us prove that 
it is true at rank t. 

1) Const (I*) C H t : Let I £ Const (I*), and let I p be the 
set of previously-detected parent targets from /, i.e., the set 

I p = {(t',y) £r t _ 1 \(f,y-z) & I,z £ Z t }. (112) 

By construction I p Q i?_1- and if | Ip | < 1 we have 
immediately I p £ Const(IJ_ 1 ). Let us suppose that \I p \ > 2, 
then let (t' ,y’), (t",y") £ I p with (t’,y') f 

where y' = (z' 0 ,..., z' t _f) and y" = (z({,... Since 

I £ Const (I*), there exist z [, z" £ Z t such that 

{t',y' :z'),(t",y" :z") £ I* and y’:z' ~ y":z". Then by 
definition of the track compatibility (65): 


Pl :z (f) = 


ao7a) 

P t \t- 

P\\t-i(9t{-, z)l m f) +pl\ t _ 1 (g t (-, z)lxj) 
Pi\t-i(9t(-,z)l w ) +p^ t _ 1 (g t (-,z) l Xf ) 

(107b) 

Pl\t-i(9t(-, z)l w f) + z)lxj) 

(1 -Pt| t -i( x t))5t(^i^) +Pt|t_i(fft(-»^)lx t )’ 

(107c) 


where we substitute Eq. (106) into Eq. (107a) to yield the 
result (107b). The probability of presence of the updated track 
i:z is then given by selecting / = lx t in Eq. (107c), i.e.. 


[z' t ,„ = z”,„ , 0 < t'" <t\^ z' t ,„ (113) 

which implies that y' ~ y". Thus, using the definition of 
consistent populations (66), I p £ Const(I*_ 1 ). Using the case 
at rank t— 1, we conclude that I p £ H t _i. Now, let us define: 

• Id = {( f,y') G Ip | C t',y':z ) £ I,Z £ Z t }- 

• Z d = {z £ Z t \ (t 1 , y' :z) £ /, (*', y') £ J d }; 

• v is the bijective function from I,\ to Z,\ such that 
(t',y':v{i)) £ I, (; t',y') £ J d ; 

• Zf d t ' = {z £ Z t | (t 1 , 4> t ~i :z ) £ I }, 0 < t r < f; 

• n t > = |Zf d ,i'|, 0 < t’ < t. 

Then, since I p £ H t _i and considering the data association 
h = (I d ,Z d ,Z fd ,n) £ Admz t (/p,n), (114) 


rf :z (x t ) = 


Pl\t-l(9t(-, Z)lyi t ) 




(108) 

Suppose that z £ Z t . Then, exploiting the stochastic descrip¬ 
tion of the sensor system (53) gives gt(f>, z) = 0, and therefore 
the probability of presence (108) becomes 


Pt z (*t) = I- 


(109) 


Suppose that z = <j>. Then, exploiting the stochastic description 
of the sensor system (53) gives gt('tj). z) = 1, and therefore the 
probability of presence (108) becomes 


P l t Z (Xt) = 


Pt|t-i(St(-)0)lx t ) 


1 -Pt| t _i(X t ) +pj |t _ 1 (5t(-,^)lx t )' 


(HO) 


then the association scheme a = (J p ,n, h) leads to the 
construction of hypothesis /; in other words, I £ H f . 

2) H t C Const(I*): Let H £ H t , and let a = (7T p ,n, h), 
where ( H p ,n ) £ H t _i x N t , h £ Adm^ t (H p , n), be the 
association scheme that produced hypothesis H (see Eq. (62)), 
with h= (//,], Z,\. Zf d . u). By construction H C I*, and if 
\H\ < 1 we have immediately H £ Const (I*). Let us 
suppose that |iT| > 2, and let (t',y'),{t",y") £ H with 

it'.y') * (t',y"). 

1. Assume y' = <p t _ 1 :z / and y' = <f t _ x \z". If z' = z" 
then y = y', and also z" £ Zf djt r and thus t" = t', which 
contradicts the assumption (t',y r ) f (t.",y"). Thus z' f z ", 
and therefore y' ~ y". 

2. Assume y' = <p t _ 1 :z' and y" = y p \z", where 
(t", I/") G H p . Then either z" = f or z" G Z d , in both options 
z’ f z" since z’ £ Zf di f, and therefore y' ~ y". The same 
reasoning applies if y' = y' p :z' and y" = 4> t _ 1 :z". 
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3. Assume y' = y' p : z' and y" = y p : z", with 

e H P- Let us P rove that (t'ly'p) ± 
for that, let us suppose that ( t',y ’) = If 

{t',y' p ) &H d , then z' = v((t',y' p )) = v((t",y'f) = z"; if 
(t ! . y' p ) £ H p \ H d , then z = <$> = z'\ in both options y' = y", 
and since ( t\y' p ) = (t” ,y”) it follows that (■ t',y ') = (■ t",y"), 
which contradicts the assumption {t',y') ^ ( t",y"). We thus 
have (t', y'p) ^ (f", y"). Using the case at rank t — 1 we have 
also {(f, y'p), (t", y'p)} Q Hp £ H t _i = Const(F t *_ 1 ) and 
the definition of the track compatibility (65) therefore yields 

<4-2/p- (H5) 

Let us now suppose that z' = z". If z' ^ </>, then 
(i', 2 /p) = ^ _1 (z') = z/ _1 (z") = (i", y"), which contradicts 
the fact that (t',y' p ) ^ (t". y"). Therefore 

[ 0 ' = 0 "] =$ z' = (f>, (116) 

and from Eqs (115) and (116) we conclude that y ~ y'. 

In all three options above we have y ~ y\ and therefore 
H £ Const (I*). 

3) H 4 = Const (I*); Since Const (I*) C H t and 
Hj C Const (I*) it follows that Ht = Const (I*), which 
proves the case at rank t. ■ 

References 

[1] Y. Bar-Shalom, “ Tracking methods in a multitarget environment,” 
Automatic Control, IEEE Transactions on, vol. 23, no. 4, pp. 618 — 
626, Aug. 1978. 

[2] D. Reid, “An Algorithm for Tracking Multiple Targets,” Automatic 
Control, IEEE Transactions on, vol. 24, no. 6, pp. 843-854, 1979. 

[3] S. S. Blackman, “Multiple Hypothesis Tracking for Multiple Target 
Tracking,” Aerospace and Electronic Systems Magazine, IEEE, pp. 5-18, 
Jan. 2004. 

[4] W. D. Blair and Y. Bar-Shalom, Eds., Multitarget-Multisensor Tracking: 
Applications and Advances (Volume III). Artech House, 2000. 

[5] G. Stuart, R. Odom, and F. Gorecki, “Implementation of a JPDAF tracker 
for an electronically scanned radar,” in AIAA Guidance, Navigation and 
Control Conference, 1989, pp. 1306-1312. 

[6] R D. Bums and W. D. Blair, “Optimal Phased Array Radar Beam 
Pointing For MTT,” in Aerospace Conference, IEEE Proceedings on, 
vol. 3, Mar. 2004. 

[7] S. S. Blackman and R. Popoli, Design and analysis of modern tracking 
systems. Artech House, 1999. 

[8] R. P. S. Mahler, Statistical Multisource-Multitarget Information Fusion. 
Artech House, 2007. 

[9] -, “Multitarget Bayes Filtering via First-Order Multitarget Mo¬ 

ments,” Aerospace and Electronic Systems, IEEE Transactions on, 
vol. 39, no. 4, pp. 1152-1178, 2003. 

[10] B.-T. Vo and B.-N. Vo, “Labeled Random Finite Sets and Multi-Object 
Conjuguate Priors,” Signal Processing, IEEE Transactions on, vol. 61, 
no. 13, pp. 3460-3475, 2013. 

[11] R. P S. Mahler, “Objective Functions for Bayesian Control-Theoretic 
Sensor Management, I: Multitarget First-Moment Approximation,” in 
Aerospace Conference, IEEE Proceedings on, vol. 4, Mar. 2003, pp. 
1905-1923. 

[12] -, “Unified Sensor Management Using CPHD Filters,” in Informa¬ 

tion Fusion, Proceedings of the 10th International Conference on, Jul. 
2007, pp. 1-7. 

[13] B. Ristic, B.-N. Vo, and D. E. Clark, “A note on the Reward Function for 
PHD Filters with Sensor Control,” Aerospace and Electronic Systems, 
IEEE Transactions on, vol. 47, no. 2, pp. 1521-1529, Apr. 2011. 

[14] M. Andrecki, E. D. Delande, J. Houssineau, and D. E. Clark, “Sensor 
management with regional statistics for the PHD filter,” in Sensor Signal 
Processing for Defence (SSPD), Sep. 2015. 

[15] M. Beard, B.-T. Vo, B.-N. Vo, and S. Arulampalam, “Sensor Control for 
Multi-target Tracking using Cauchy-Schwarz Divergence,” in Informa¬ 
tion Fusion, Proceedings of the 18th International Conference on, Jul. 
2015, pp. 937 - 944. 


[16] J. Houssineau, “Representation and estimation of stochastic popula¬ 
tions,” Ph.D. dissertation, Heriot Watt University, 2015. 

[17] E. D. Delande, M. Uney, J. Houssineau, and D. E. Clark, “Regional Vari¬ 
ance for Multi-Object Filtering,” Signal Processing, IEEE Transactions 
on, vol. 62, no. 13, pp. 3415-3428, 2014. 

[18] E. D. Delande, J. Houssineau, and D. E. Clark, “Performance metric in 
closed-loop sensor management for stochastic populations,” in Sensor 
Signal Processing for Defence (SSPD), Sep. 2014. 

[19] R. E. Kalman, “A New Approach to Linear Filtering and Prediction 
Problems,” Journal of Basic Engineering, vol. 82, no. 1, pp. 32—45, 
1960. 

[20] P. Del Moral and A. Doucet, “Particle methods: An introduction with 
applications,” INRIA, Tech. Rep., 2009. 

[21] D. Stoyan, W. S. Kendall, and J. Mecke, Stochastic geometry and its 
applications, 2nd ed. Wiley, 1995. 

[22] D. J. Daley and D. Vere-Jones, An Introduction to the Theory of Point 
Processes, 2nd ed., ser. Statistical Theory and Methods, D. J. Daley and 
D. Vere-Jones, Eds. Springer Series in Statistics, 2003, vol. 1. 

[23] -, An Introduction to the Theory of Point Processes, 2nd ed., ser. 

Statistical Theory and Methods, D. J. Daley and D. Vere-Jones, Eds. 
Springer Series in Statistics, 2008, vol. 2. 

[24] J. Houssineau, D. E. Clark, and P Del Moral, “A sequential Monte 
Carlo approximation of the HISP filter,” in 2015 23rd European Signal 
Processing Conference (EUSIPCO), Aug. 2015, pp. 1251-1255. 

[25] D. L. Alspach and H. W. Sorenson, “Nonlinear Bayesian estimation 
using Gaussian sum approximations,” Automatic Control, IEEE Trans¬ 
actions on, vol. 17, no. 4, pp. 439^448, 1972. 

[26] A. Doucet, N. de Freitas, and N. Gordon, Sequential Monte Carlo 
Methods in Practice, ser. Statistics for Engineering and Information 
Science. Springer, 2001. 

[27] A. Renyi, “On measures of entropy and information,” in Fourth Berkeley 
Symposium on Mathematical Statistics and Probability, vol. 1, 1961. 

[28] V. I. Bogachev, Measure Theory. Springer, 2007, vol. 1. 

[29] E. D. Delande, J. Houssineau, and D. E. Clark, “A filter for distinguish¬ 
able and independent populations,” 2015, arXiv: 1501.0467lvl. 

[30] E. D. Delande, C. Frueh, J. Franco, J. Houssineau, and D. E. Clark, “A 
Novel Multi-Object Filteimg Approach For Space Situational Aware¬ 
ness,” Journal of Guidance, Control, and Dynamics, 2016, submitted. 


